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Abstract 

We study the critical behavior of the 3-state Potts model, where the spins are 
located at the centers of the occupied squares of the deterministic Sierpinski car- 
pet. A finite-size scaling analysis is performed from Monte Carlo simulations, for 
a Hausdorff dimension df — 1.8928. The phase transition is shown to be a second 
order one. The maxima of the susceptibility of the order parameter follow a power 
law in a very reliable way, which enables us to calculate the ratio of the exponents 
j/u. We find that the scaling corrections affect the behavior of most of the thermo- 
dynamical quantities. However, the sequence of intersection points extracted from 
the Binder's cumulant provides bounds for the critical temperature. We are able to 
give the bounds for the exponent 1/V as well as for the ratio of the exponents fijv, 
which are compatible with the results calculated from the hyperscaling relation. 

Keywords: phase transitions, fractal, Potts model, finite size scaling, Monte Carlo. 
PACS: 68.35.Rh: Phase transitions and critical phenomena, 05.45.Df: Fractals, 75.10.Hk: 
Classical spin models, 75.40.Mg: Numerical simulation studies, 89. 75. Da: Systems obey- 
ing scaling laws. 

1 Introduction 

Critical phenomena on fractal structures have been first studied by Gefen et al. pQ. 
Since then, most of the works on this topic have been done in the framework of the Ising 
model. Recently, the use of cluster algorithms and histogram methods enabled two groups 
[3 El lU El to study the ferromagnetic Ising transition on the fractal lattices by Monte 
Carlo simulation much more thoroughly than before; several fractal dimensions between 
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1 and 3 have been investigated. These last results show that the scaling corrections can 
strongly affect the finite size behavior of some thermodynamical quantities, especially 
when the fractal dimension is lowered from 2 to 1; the critical exponents cannot always 
be calculated, but bounds can be provided. The available values and bounds of the 
exponents are consistent with the hyperscaling relation provided that the space dimension 
is replaced by the Hausdorff one. Furthermore, they do not agree with e-expansion results 
jSj. The critical behavior of fractals is said to be understood in the framework of weak 
universality [UHllZj: the critical exponents do not only depend upon the symmetry of 
the order parameter, the interaction range, and the fractal dimension, but also upon 
geometrical features of the fractal structure. 

The Potts model is a generalization of the Ising model obtained by varying the number 
of spin states from 2 to any non zero value q. Since the number of states q is related to the 
symmetrical properties of the order parameter, the critical behavior of the Potts model 
depends upon the value of this additional variable One of the most striking results is 
the effect of q on the order of the transition as a function of the space dimension d: In the 
case of regular translational invariant lattices, it has been shown that the ferromagnetic 
phase transition is a second order one if q is smaller than a value q c (d) and a first order 
one if q > q c {d). For instance, the phase transition in the bidimensional case is first order 
for q > 4 and second order for q < 4; the value of q c (3) is smaller than 3. Moreover, 
a wide variety of effects occur in the presence of disorder: bond randomness can induce 
a second order transition from a system exhibiting a first order one 0- The order of 
the transition can be changed to a second one if strong enough aperiodic fluctuations are 
introduced [TUj . 

The question arises whether the values q c (d) make sense in non integer dimensions dj 
or not. Fractals are natural candidates to interpolate between integer space dimensions. 
Since translation invariance is broken in these structures, disorder is introduced: the 
lattice is strongly inhomogeneous. Thus, two among the three parameters driving the 
critical behavior of the Potts model (structured fractal disorder and dimensionality) are 
closely linked in the present case of fractals. The critical properties of the Potts model on 
fractal structures have never been studied before, excepted by Bin [TTj who generalizes 
the Migdal-Kadanoff bond-moving renormalization scheme on the Sierpinski carpet from 
the Ising model |T2] to the Potts one. No Monte Carlo results are available, up to now; 
there are two main reasons why studies of the critical behaviors on fractal structures is a 
difficult task: 

(i) Monte Carlo simulations come up against critical slowing down. Fortunately, the 
use of cluster algorithms [T3*l ITi] is very helpful in overcoming this difficulty 

(ii) The scaling corrections are expected to be large [5]. A reliable analysis, hence, 
requires simulations on several large lattices of the same structure. Due to the 
self-similar character of the underlying network, the sizes of the lattices increase 
as a geometrical series. We face the high C.P.U. time consuming problem in the 
simulations. 
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The purpose of this paper is to provide a thorough study of the 3-state Potts model 
on the deterministic Sierpinski carpets of a given fractal dimension. The order of the 
transition has to be carefully checked before calculating the critical temperature and 
exponents. The paper is organized as follows: The model and the finite size scaling 
theory are briefly recalled in Sec. 2.1 and Sec. 2. 2. The simulation methods are described 
in Sec. 2. 3. The critical behavior of the 3-state Potts model on the deterministic Sierpinski 
carpet is studied in Sec. 3. 

2 Methods and theoretical background 

2.1 Sierpinski carpet SC(3,1) and Potts model 

The lattice structures we deal with are constructed according to an iterative segmen- 
tation process: we start from a generating cell denoted SC(3, 1,1), which is a 3 x 3 square 
lattice where the central square is removed. The lattice SC(3, l,k) associated with the 
k th iteration step is constructed by replacing each occupied square of SC(3, 1, k — 1) by 
the generating cell while enlarging the side 3 times. The size of the lattice SC(3, 1, k) is 
L = 3 k and the number of occupied sites iV£ cc = (3 2 — l 2 ) k . The Hausdorff dimension 
df, defined in such a manner that iV£ cc = L f, is equal to ln(8)/ln(3) m 1.892789. The 
Sierpinski carpet SC(3,l,k) becomes a true fractal in the mathematical sense when k 
tends to infinity, and we will denote it SC(3, 1). 

The Hamiltonian of the g-state Potts model on SC(3, 1, k) reads : 



where <jj and <jj stand for the spin states on the sites i and j and can take the integer 
values 1,2, ... ,q. 5{ai,aj) is the Kronecker ^-function, and J is the coupling constant. 
The sum runs over the nearest-neighbor bonds on SC(3, l,k). The order parameter per 
site of the Potts model is defined as : 



where pi = max(ni, ri2, . . . , n q )/N^ cc is the largest density of spin state, since n s is the 
number of spins in state s. We notice that the values of lie in a range between and 
1. We will focus our attention on the 3-state Potts model. 

2.2 Finite size scaling for second order and first order phase 
transitions 

Finite size scaling theory has been first developed by Fisher et al. [TS] in the case 
of the second order phase transitions. According to Widom's homogeneity hypothesis, 
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applying to a system with translational symmetry, the singular part of the free energy 
per spin of the system under a change of the unit length from 1 to b is assumed to scale 
as 

f(t,h,L) = b- d f(W,hby\L/b) (3) 

where d is the lattice dimension, L is the original lattice size, yt and yn are eigenvalue 
exponents, t = (T — T c )/T c is the reduced temperature with T c representing the critical 
temperature on the infinite lattice, and h is the external field; t and h are supposed to 
be small, and b smaller than the correlation length £. This hypothesis (Eq.(J5|)) has to 
be generalized in the case of fractals by restricting the choice of b to a power of the size 
of the generating cell of the Sierpinski lattice in order to keep invariant the structure 
of the renormalized lattice. Also, one has to replace the dimension d in Eq.(j3J) by the 
Hausdorff dimension df of the fractal; the factor b f describes how the number of spins 
decreases under the change of the unit length. Moreover, since the translational symmetry 
is lost on Sierpinski lattices, the two point correlation function depends on the positions 
of the spins. Hence, the definition of the correlation length on the fractal lattices has 
to be adapted. A position-independent correlation function of one scalar variable r can 
be defined by taking the average of the two point correlation function values over all 
the pairs of spins having direct distance equal to r. The position- independent correlation 
function behaves as exp(— r/£(T))/r d_2+7? in the critical region in the case of translational 
symmetry networks. We assume that it keeps a similar behavior on fractal lattices too. 
For a second order phase transition system, £(T) diverges at the critical point. We, 
therefore, expect that £(T) reaches the finite lattice size L while T is in the neighborhood 
of T c . In this temperature region, we can take b equal to the limiting value L and it yields 
the zero-external field scaling equations for the fractal lattice as: 

C L (T) ~ L a/v TC(tL 1,l/ ) (4) 

m L (T) ~ L~ p/V M(tL 1/u ) (5) 

X l(T) ~ L^X(tL^) (6) 

L (T) ~ L l/u T 2 V(tL 1/v ) (7) 

U L (T) ~ TU(tL 1/u ) (8) 

where Ci(T) is the heat capacity per site, mz,(T) the order parameter per site, Xl(T) the 
zero field susceptibility of the order parameter per site, 4>l(T) = 91n(mz,(T)) /d(ksT)^ 1 
usually called the first logarithmic derivative of the order parameter, and Ul{T) = 1 — 
(m| / (T))/(3(m| / (T)) 2 ) the Binder's cumulant. The critical exponents a, (3, 7 and v can 
be written as a function of the two independent eigenvalue exponents y t and y^ and 
the dimension df. a is equal to 2 — df/y t , [3 is equal to (df — yh)/yt, 7 is equal to 
(2yh — df)/y t , and v is equal to 1/yt- The finite- size shifting of the rounded singularity 
of a given thermodynamical average k (for instance, 0, Xi an d C) from the critical point 
should follow the relation : 

T%L)=T c + g K L-^, (9) 
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where the g K 's are some constants. Moreover, the width of the associated rounded singu- 



For a system undergoing a first order phase transition, Fisher et al. JB] and Cardy 
et al. [T7] propose that the relevant eigenvalue exponents in Eq.® should reach their 
allowed limiting values, said y t = d and y^ = d, in the phenomenological studies. In the 
case of the Sierpinski lattice, yt and yh should take the value of the Hausdorff dimension 
df. The thermodynamical averages on the fractal lattice, therefore, scale as: 



Moreover, they are characterized by 5-function like singularities as the lattice size L 
increases ; as a result, the areas under the curves displaying a given average as a function 
of the temperature should be independent of L. Also, it is well known that the probability 
distribution of the energy is double peaked in the vicinity of the transition temperature 



2.3 Simulation process 

We implement the Wolff 13 j and Swendsen-Wang Monte Carlo cluster algo- 
rithms; both of them have been shown to largely overcome the critical slowing down. In 
order to save simulation time, we used alternatively the two algorithms in our studies. 
Let us describe the simulation process in a more precise way. An initial configuration is 
built up by randomly assigning one of the three spin states to each occupied site of the 
lattice SC(3, l,k). Periodic boundary conditions will be used. A simulation temperature 
To is set and enough Wolff or Swendsen-Wang Monte Carlo steps are performed in order 
to ensure that the configuration has been thermalized. The energy per site eL(To) and 
the order parameter per site m L (T ) are then calculated for every Monte Carlo step. One 
million steps are collected to be a sample and the above procedure is started again at 
least ten times. The histogram method j2Uj . applied within each sample, allows the ther- 
modynamical averages of a physical quantity g to be calculated over a range AT around 
T ; we denote it (g(T ))T- The reliability of the calculated averages over the temperature 
interval must be carefully checked on. <p L (T), xl(T), Cx(T), and Ul{T) are then obtained 
according to the following relations: 



larity should scale as L x l v . 



Cl(T) 
Xl(T) 



L d fTC{tL d f) 
L d f X(tL d > ) 
L d f T 2 V(tL d f) 



(10) 

(11) 
(12) 
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Xl(T) 



N°~((ml(T )) T - (m L (T ))l) 
k B T 



(14) 
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Ul(T) 



1 - 
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3(m|(T ))|- 



(16) 



The statistical errors associated with the thermodynamical averages are estimated by 
calculating the standard deviations from the ten samples. 

3 Results and discussions 

3.1 Studies of rounded singularities of Xl(T) and 4>l(T) 

First, we study the order of the phase transition of the 3-state Potts model in the 
case of the fractal SC(3, 1). The behavior of the finite-size rounded singularities of Xl{T) 
calculated from our simulations are shown in Fig.l for L = 27, 81, 243, 729, and 2187. 
The maximum values xT ax an d their associated temperatures T X (L) are reported in Table 
1. The finite-size shifts of XL ax are clearly observed; the Lorentz distribution is a better 
candidate to fit the rounded singularity than the Gaussian one: 



The width of the Lorentz distribution is W and the area under the distribution curve is 
A. The two-parameter fits (W, A) of the rounded singularities are reported in Table 2. It 
appears clearly that the area A under the Lorentz distribution increases as L increases; 
as a result of this, we can conclude that the singularity is not 5-function like as L tends 
to infinity. We have the first evidence that the phase transition is not a first order one. 
According to the finite size scaling for second order phase transitions, the width W should 
scale as L~~ l l u and the area A as L 1 ^^ 1 ^ . The log-log plots of W and A versus L do 
not appear straight lines over the whole range from L = 27 to L = 2187; they exhibit a 
slight curvature, which can be interpreted as a contribution of scaling corrections to the 
behavior of W and A. The slopes calculated from two consecutive points from L = 27 to 
L = 2187 are -0.786(17), -0.701(11), -0.6243(60), -0.5524(35) for W and 0.897(12), 
0.9867(81), 1.0635(40), 1.1634(25) for A. The effect of the scaling corrections in the 
case of fractal lattices has been pointed out and discussed in the framework of the Ising 
model by Monceau et al. [5 j and Carmona et al. jHJ. Instead of a power law, W is 
expected to behave as L~ x l v {\ +(pw(L)) where the scaling corrections <p w (L) tend to zero 
as L tends to infinity; usually, <pw(L) is developed as (pw(L) = AwL~ u , where Aw is 
the amplitude of the corrections and u an additional exponent. A similar behavior is 
expected for A. There is no hope to calculate the scaling corrections in a reliable way 
from our simulations. However, the slopes calculated above enable to estimate an upper 
bound for \jv to be 0.5524(35) and a lower bound for 7/1/ — 1/u to be 1.1634(25). 

We now pay attention to the maxima Xi ax °f Xl(T), and 0™ ax of <Pl(T), reported in 
Table 1. The log-log plot of XL ax shows a straight line. The power law fits are satisfied 
with reliability coefficients equal to 1.0000; they yield 7/1/ = 1.6948(37) with a four sizes 



Xl(T) 



2A r w 

~ [A(T -Tx(L)) 2 + w 2 



(17) 
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fit from L = 81 to L = 2187 and 7/z/ = 1.7013(28) with a three sizes fit from L = 243 to 
L = 2187. The scaling corrections turn out to be very weak in the case of the maxima of 
the susceptibility, as already shown in the case of the Ising model. The best estimate of 
7/z/ is 1.7013(28). Finite size scaling theory states that the slope of the line represents the 
exponent 7/1/ if the phase transition is second order and that it should be the Hausdorff 
dimension df if the transition is first order. The stability of the fits involving xT ax 
enable to conclude that 7/1/ is significantly different from df = 1.892789 and, therefore, 
consistent with a second order phase transition as we have shown at the beginning of 
this paragraph. On the other hand, 0™ ax suffers from large scaling corrections since 
the slopes of two successive points of <p™ ax in a log-log plot are measured as 0.758(20), 
0.669(11), 0.613(13), and 0.533(11). These large scaling corrections may be due to the 
deviation of the average number of nearest neighbors on the finite size fractal from that 
on the infinite fractal [5] because the critical temperature is directly related to the average 
number of nearest neighbors in the mean filed theory and the eignvalue exponent y t = 1/v 
describes the scaling behavior of the reduced temperature. The behavior of 0™ aa: provides 
an upper bound for 1/V, which is 1/v — y t < 0.533(11). Combining 7/z/ = 1.7013(28) 
with 1/v < 0.533(11), we get a lower bound for (7/z/ — 1/v), which is 1.168(11). These 
bounds are both very close to the ones we got from the fits of W and A. The eigenvalue 
exponent can be calculated from (7/z/ + df)/2 and is equal to 1.7970(15). 

The comparison between the present value of 7/z/ and the one associated with the 
Ising model for the same fractal structure (7/z/ = 1.732(4)) |3] enables to distinguish the 
Ising and the 3-state Potts model universality classes. We mention that this comparison 
makes sense only if the fractal structures are exactly the same because the universality 
is said to be weak [TJ HJ [7j . 



3.2 Estimation of the critical temperature T c 

The infinite lattice critical temperature T c can be calculated from the finite-size shifts 
of the rounded singularities provided that v is known; for a given peaked thermodynamical 
average k, the equation Q yields: 

T"{lL)-rV»T"{L) 
Tc = 1 _ £-1/, ( 18 ) 

where i = 3 is the size of the generating cell. In this form, the value of T c is deduced 
from two consecutive iteration steps of the fractal structure. An upper bound for 1/v 
will yield an upper bound for T c since T X (L) and T^(L) decrease as L increases. Taking 
1/v = 0.533(11), the values of T X (L) from the two highest iteration steps provide the 
first estimate for the upper bound of T c equal to 0.67324(1264). The values of T^(L) 
from the two highest iteration steps yield 0.67357(1456) as a second upper bound for T c . 
The errors in the above bounds mainly come from the error propagation of 1/v. We can 
conclude that T c < 0.673(12). 

An alternative way to measure T c without knowing v is to study the Binder's cumulant 
crossing points |19] . Fig.2 shows that the cumulant curves for different lattice sizes do 
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not intersect at a fixed point. We call Tfc-i^ the temperature of the cumulant intersection 
point between the iteration steps (k — 1) and k. The measured values are: T 2 3 = 
0.68461(54), T 3i4 = 0.67923(26), T 4>5 = 0.67594(20), T 5j6 = 0.673916(69), and T 6J = 
0.672866(50). The difference AT^ = (T/ C _ lj fc — fc+i) behaves basically as a geometrical 
series as a function of k. If we write AT^ = ar 3 r 4 ■ • ■ where = AT^/AT^!, k > A, 
and ar 3 = AT 3 , we will find 1 > r 4 > r 5 > r 6 > • • •. This enables us to claim that 
the series {T^+i} converges as k tends to infinity. The convergent value calculated 
from the three intersection points T 4i 5, T^^, and Tqj gives a lower bound for the critical 
temperature: 0.671737. We finally conclude that 0.671737 <T C < T 6J = 0.672866. This 
range is in agreement with the upper bound for T c extracted from Eq.(|9*jl. 

3.3 Simulation near T c 

We firstly check the occurence of a single peak in the energy probability distributions, 
which confirms again that the transition is second order. A typical energy probability 
distribution is shown in Fig. 3. We then study the behaviors of 4>l(T), Xl(T), and m^T) 
within the previously estimated temperature range for T c . The following three tempera- 
tures will be considered: 0.671737, 0.672866, and 0.673241. The first two temperatures 
are the lower and the upper bounds extracted from the Binder's cumulant crossing and 
the last one is the upper bound calculated from Eq.([18|). assuming that \jv — 0.533 and 
using the positions of the peaks of xl at L = 729 and L = 2187. The log- log plots of 
mi(T), xl(T), and </>l(T) at the three considered temperatures are shown in Fig.4. mz, 
shows a straight line in the log-log plot, at T = 0.671737, over a range of lattice sizes 
covering several orders of magnitude, with a slope, measured from L = 243, L = 729, 
and L = 2187, equal to —0.0748(1). When the temperature increases towards the upper 
bound, mi(T) becomes more and more concave while Xl{T) is more and more convex. 
The asymptotic slopes for m^(T) obtained from the measurements beween the largest 
sizes L = 729 and L = 2187 are upper bounded at T = 0.672866 and T = 0.673241, 
respectively, by —0.0876(2) and —0.0932(2); and the asymptotic slopes for Xl(T) at 
T = 0.671737, T = 0.672866, and T = 0.673241 are lower bounded by 1.60(1), 1.73(1), 
and 1.78(1). On the other hand, the curvature of <Pl{T) in a log- log plot varies in the 
opposite direction, that is, the curve is less convex at the higher temperature bound 
(T = 0.673241) than at the lower one (T = 0.671737). The asymptotic slopes for 4>l(T) 
at T = 0.671737, T = 0.672866, and T = 0.673241 are upper bounded by 0.466(13), 
0.521(10), 0.552(11). The scaling corrections, depending on the physical quantities, turn 
out to affect strongly the critical behavior of the 3-state Potts model on the fractal struc- 
ture even when the sizes are large. The estimation of the ratio of critical exponents [3/v, 
7/V, and \jv from the behavior of rriL, Xl, and <pL in the critical region becomes a very 
difficult task: firstly, the value of the critical temperature cannot be provided precisely; 
secondly, even if we expect scaling corrections to tend to zero as L tends towards infinity 
(whatever thermodynamic quantity is studied), they should be included in the calcula- 
tion of our exponents; as already explained in Sec. 3.1, there is no hope to extract scaling 
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corrections from our data. We point out that the scaling corrections in the critical region 
are stronger than that for the case of Ising model. In the case of Ising model, the slope of 
Xl{T) extracted at the expected critical temperature is consistent with the value of 7/1/ 
extracted from x™ ax El E] ■ 

3.4 Discussions 

The consistency of our full set of results can be achieved in the following way: The 
behavior of the Binder cumulant shows that T c should be greater than 0.671737 and 
smaller than 0.682866. From the behavior of xl{T), the asymptotic slopes in this tem- 
perature region are lower bounded by 1.60(1), consistent with the value j/u — 1.7013(29) 
extracted from x r E ax - The asymptotic slopes for m^(T) in the region are always smaller 
than -0.0748(1) and, thus, yield (3/v > 0.0748(1), consistent with the requirement of the 
hyperscaling relation, said /3/u = (df — r y/v)/2 = 0.0957(21) with 7/V taking the value 
1.7013(29). Similarly, 1/v < 0.521(11) can be found from the behavior of 4>l{T) in this 
temperature region, consistent with the results in the studies of (f)™ ax . We, therefore, give 
the best value and bounds of the exponents and the critical temperature that we find: 
7/1/ = 1.7013(29), (3/v > 0.0748(1), 1/u < 0.521(11), and 0.671737 <T C < 0.672866. 

The critical temperature and the exponents of the 3-state Potts model on a two 
dimensional regular system are known to take the following values |2I]: T c = l/ln(l + 
v 7 ^) ^ 0.995, l/v = 6/5, p/u = 2/15, and j/p = 26/15. We find that the values of T c , 
1/v, (3/v and 7/1/ for the 3-state Potts model are smaller in the case of the Sierpinski 
carpet we investigated than in the case of the two dimensional regular system. A similar 
situation occurs in the case of the Ising model on Sierpinski carpets .5] where, regardless of 
the structure of the fractals, those values decrease as the Hausdorff dimension decreases. 

Discrepancies are observed between our simulation results and Bin's predictions. Bin 
[TT] found by applying the method of Migdal-Kadanoff bond-moving renormalization, 
that exp(— K) = 0.746 (which is equivalent to T c = 3.4126) and i/k(— Vt) = 0.619 for 
the 3-state Potts model on the same Sierpinski carpet we studied in this paper. There 
are mainly two reasons for the occurence of discrepancies: firstly, the Migdal-Kadanoff 
bond-moving renormalization is an approximate method; secondly, in Bin's case the spins 
were located at the vertices of the occupied squares of the Sierpinski carpet while they 
are located at the centers of the occupied squares in the present work; it yields a different 
universality class because the duality is no more held on fractal lattices. 

In conclusion, we have carefully performed the first Monte Carlo studies for the 3-state 
Potts model on the Sierpinski carpet with df ~ 1.8928. Although the system exhibits 
strong scaling corrections, the consistency of our results demonstrates the validity of the 
finite size scaling theory for the g-state Potts model on these fractal lattices. Nevertheless, 
scaling corrections tend to increase while q is passing from 2 (Ising model) to 3. We, hence, 
expect the raise of these difficulties in the studies of the phase transition while q goes 
larger. In the coming future, we will be able to perform in a reliable way the studies for 
the Potts model with larger q on fractal dimensions higher than 2. 
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Figure captions 

Fig.l Rounded singularities of the susceptibility on lattice SC(3, 1, k) where k — 3, 4, 5, 
6, and 7. Each segment represents the reliabe temperature range of the histogram 
method where the simulation temperature is located in its center. The fitting by 
Lorentz distribution for the susceptibility on 5(7(3, 1, 7) is plot in order to clearly 
indicate the location of the rounded singularity. 

Fig. 2 Binder's cumulant crossings. 

Fig. 3 Probability distribution of energy per site on 5(7(3,1,7). It's made of 10 million 
Monte Carlo data and seperated by 500 histogram classes. The simulation temper- 
ature is taken at 0.672845. 

Fig. 4 Log-log plot of rriL, Xl, and 0^ at T = 0.671737 (open square), 0.672866 (open 
circle), and 0.673241 (open triangle) on the Sierpinski carpet 5(7(3, l,k). 
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Table 1: Maxima values and associated peak positions of xl and 0l on the Sierpinski 



carpet 5(7(3, 1, k) of size L = 3 k 





L = 27 


L = 81 


L = 243 


L = 729 


L = 2187 


max 

Xl 

TX(L) 


30.53(12) 
0.73800(22) 


195.58(59) 
0.70448(12) 


1239.4(4.4) 
0.689525(65) 


7984(19) 
0.682115(16) 


51969(105) 
0.678182(11) 


jjnax 

<Pl 

T*{L) 


8.85(18) 
0.7510(17) 


20.35(13) 
0.70995(38) 


42.46(45) 
0.69221(22) 


83.23(74) 
0.68357(16) 


149.50(1.23) 
0.679139(23) 



Table 2: Fit of the Lorentz distribution (Eq. (fT7j) ) of the rounded singularity of xl on the 



Sierpinski carpet SC(3, l,k) of size L = 3 k . 





L = 27 


L = 81 


L = 243 


L = 729 


L = 2187 


w 

A 


0.0700(11) 
3.444(38) 


0.02949(31) 
9.230(74) 


0.01366(8) 
27.29(11) 


0.00688(2) 
87.78(18) 


0.00375(1) 
315.13(57) 
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